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PREFACE 

THE aim of the present volume is to set forth a method of treating 
the geometrical optics of any given system intrinsically, that is to 
say, without introducing any artificial scaffolding and, therefore, 
without any arbitrary splitting of the entities involved in the 
problem of tracing a luminous ray through the several surfaces 
of a system. The resultant formulae, in vector language, of 
course, as the only appropriate means of intrinsic expression, will 
then be free of any unnecessary elements, and therefore, if finally 
translated into any set of scalar formulae for the practical use 
of the numerical computer, will contain no superfluous, geometrical 
or arithmetical, complications. 

Our purpose is not to treat the whole subject of geometrical 
optics, but exclusively, or almost so, that part of it which is called 
by the short name of " ray tracing. 1 ' This is notoriously the most 
laborious part of the computer's patient work, and becomes, 
without question, a formidable task when he has to deal with 
skew rays and non-centred systems. The problem thus limited 
can be put shortly : Given the ray incident upon any system 
of lenses, mirrors, and prisms, find the emergent ray. 

The advantages of the vectorial method of resolving it must 
not be judged by the conspicuous shortness of the resultant formulae 
alone, but also by the simplicity of their deduction, as compared 
with the usual method, and by the facility of recalling the formulae, 
or of reconstructing them if forgotten. Again, the nature of the 
proposed method is such as to make the help of drawings which, 
especially in the case of skew rays, become, in the best standard 
treatises, exceedingly complicated almost superfluous. This 
circumstance will be particularly welcome to readers who are not 
endowed with strong visualising powers. In connection therewith 
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the often troublesome discrimination of sign will become almost 
automatic. 

As already indicated, the mathematical idiom to be employed 
to secure these advantages and to prevent the intrusion of foreign 
and artificial elements is Vector Algebra. The reader whose 
acquaintance with that natural language of space-relations is but 
slight need not, on that account, be deterred from studying the 
methods to be described. For they will involve only the rudiments 
of vector algebra : namely, the addition and the scalar multiplica- 
tion of vectors together with an occasional application of the 
vector product of two vectors. A reader who is entirely ignorant 
of the subject could, in a few hours, acquire this amount of know- 
ledge from one of the existing treatises on it. For instance, the 
first twenty pages of the author's Vectorial Mechanics (Mac- 
millan, 1913) contain all that is required and even more. It is 
unnecessary, perhaps, to point out that the rudiments of vector 
addition and multiplication, thus readily learnt, will be found 
useful in other connections, such as the study of electro-magnetism 
and kinematics. It is, however, permissible to urge that British 
students, above all others, should be more widely acquainted with 
the language of vectors, so powerful and yet comparatively so little 
used. For it is their own patrimony. After due acknowledgment 
has been paid to Grassman it remains true that the vector method 
originated in this country through the creative genius of the 
great Hamilton, and was shaped into a practical instrument of 
investigation by Oliver Heaviside's pioneer work of simplification, 
seconded in America by the independent and equally successful 
efforts of J. Willard Gibbs. 

In order to avoid digressions in the text of this volume it will 
be well to give here a few explanations of the symbols to be used 
and to collect the very meagre number of vector formulae which 
will be required for our purpose. Vectors will be printed in 
Clarendon, small or capital type. Thus 

A, r, etc. 

will be vectors. Their sizes or absolute values will be denoted 
by the corresponding Italics, as A, r in the above case. If r = l, 
then r is called a unit vector. The vector sum of two vectors 
will be denoted, as usual, by A+B. The scalar product of two 
vectors A, B will be written simply AB, and their vector product, 
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which will not often be needed, VAB. With these symbols the 
only formulae required will be : 

A+B=B+A; AB = BA; VAB = -VBA, 
A(B+C)=AB+AC; VA(B +C)=VAB+ VAC. 
Similarly for the product of two sums of any number of vectors. 
If 7 is the angle contained between the two vectors A, B, then 

AB = AB cos 7, 

and the size of the vector product VAB is 
| VAB | = AB sin 7. 

The scalar autoproduct of a vector A, or its " square," will be 
written AA = A 2 = A 2 . 

If r is a unit vector, it is evident that r 2 = I. The angle 7 between 
two unit vectors a, b, is given by 

ab = cos 7. 

The mutual perpendicularity or orthogonality of any two 
vectors A, B is equivalent to AB o, while their parallelism is 
expressed by VAB=o. Instead of the latter the form A || B will 
sometimes be employed. The concurrency (or equal sense) of 
parallels and its reverse will, whenever needed, be expressly 
stated. Finally, for any A, B, 0, 

VAVBC =B(CA) - C(AB). 

Scarcely more than what has here been set down will be required 
for following freely the course of the deductions. 

The main subject will, now and then, be illustrated by examples 
or exercises. These will, for the sake of better discrimination be- 
tween essentials and accessory matter, be printed in small type. 

I gladly take the opportunity of expressing my best thanks 
to Mr. Frank Twyman, Manager to Messrs. Adam Hilger, Ltd., 
for encouraging this little work, and to my friend Prof. T. Percy 
Nunn for reading and revising the diction of the MS. My 
thanks are also due to the Publishers for the care they have bestowed 

on the book. 

L. SILBERSTEIN. 

LONDON, April 1918. 
RESEARCH DKPT., ADAM HILGER, LTD. 
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SIMPLIFIED METHOD OF TRACING RAYS THROUGH 
ANY OPTICAL SYSTEM. 

1. The Fundamental Reflection and Refraction Laws. 

LET the unit vector n represent the local normal to the surface 
separating two homogeneous isotropic media whose refractive 
indices are p. and /A'. Whether n is drawn towards the first or 
towards the second medium, is immaterial, as will be seen presently. 
Let the direction of the incident ray be given by the unit vector r, 
that of the refracted ray by r', and that of the reflected one by r", 
the latter two being again unit vectors. The plane r, n is then 
the incidence plane. The angles of incidence i, of refraction i', 
and of reflection i" ', Fig. i, are given, according to the very definition 




FIG. i. 

of the scalar and the vector product, by 

cos i =rn, cos i' =r'n, cos i" =r"n, 

or by the sizes or absolute values of the corresponding vector 
products, viz. s in i =| Vrn | , 

and similarly for the dashed angles. Thus the law of reflection, 
which asserts the equality of i, i" and the coplanarity of r, n, r", 
can be written, regard being paid to the sense of r", 
Vr'n^Vrn, r"n=rn.* 

* The second equation, besides being obvious, is also a consequence of the first. 
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Multiply both sides of the first equation vectorially by n. Then r 
by the last formula of the Preface, and remembering that n 2 = i, 
the vector form of the law of reflection will become 

r'-r-2n(nr), (I) 

giving the reflected ray r" in terms of the incident ray r and the 
incidence normal n. 

Again, the refraction law, which asserts that /*'sin I'=/A sin i 
and that r, n, r' are coplanar, is tantamount to writing 

Vns'=Vns, 

where s and s' stand for /ir and /x'r', respectively, and might, 
therefore, be called directed light paths. The name is of small 
importance. Keep in mind, however, that the r's being unit 
vectors, B 2 =/x 2 and s' 2 =/i' 2 . Treat the last equation exactly as 
in the case of reflection, i.e. multiply it vectorially by n. Then the 
result will be 8'=8+n(s'n-sn), 

or, putting g =s'n -sn =/*' cos i' -p cos i 

as a convenient abbreviation, to be used throughout the book, 

S' = B+*H. (II) 

This is the required vector form of the law of refraction,* giving 
the refracted ray s' in terms of the incident s and the normal n. 
It is true that the scalar, or ordinary numerical, coefficient g itself 
contains the unknown vector s', but only in the combination 
s'n /A'COS i'. Now n, r being given, i is given; whence i' 9 and 
therefore also g, can be rapidly calculated by means of the relation 
ft'sin t'=ft sinz. It is, therefore, not advisable to eliminate 
s' formally from, the right-hand member of (II), although this 
can easily be done. (Vide infra, Section on Multiple Reflection.) 
It will be found more convenient to leave the form (II) of the 
refraction law with the numerical coefficient g as it stands. 

* Formula (II) has also been arrived at, in a roundabout way, by Sommerfeld 
and Runge, Ann. d. Physik, 35, 1911, p. 277, the object of the authors being to 
give a general differential basis to geometrical optics. In our case it appears 
simply as the immediate vector translation of the law of refraction considered as a 
given law. Cartesian (three formulae) equivalents of (II) are, of course, to be 
found also in some of the old works. Formula (I) and its consequences were 
given, in dyadic form, in my paper "On Multiple Reflexion," Phil. Mag. for 
November 1916. 
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A glance at (I) will suffice to see that the law of reflection follows 
from that of refraction, (II), by putting in g 

-p. instead of /*', (III) 

a well-known rule. In fact, it is enough for this purpose to multiply 
(I) on both sides by the refractive index /i of the medium common 
to r and r". Remembering the simple rule (III), we can henceforth 
use (II) alone as the more general formula. Whenever a separating 
surface is contemplated as a reflector write g= -2sn = -2/* cos i. 
Unless the contrary is expressly stated we shall conceive every 
such surface as a refracting one. 

Notice that if in (II), and therefore also in (I), n is replaced by 
-n, the formulae remain intact. . 

Let now the refracted ray impinge upon another refracting 
surface. Then the new refracted ray will be obtained from the 
ray s' as this was from s, and so on, for any number of surfaces. 
Call n lf n 2 , etc., in general n*, the unit incidence normals of the 
consecutive surfaces (at the places of incidence of the ray in 
question) and /**, /*/ the indices of the media separated by the 
*th surface (Fig. 2), so that ^'=/ 



etc. 



FIG. 2. 

Further write, for each surface, 

B = J*A, a/ = /^'ir/, 

as above for a single surface, and put g K =^B K f n K -B K n K . Then, 
applying the law (II) to the successive surfaces and remembering 
that S^SK+I (newly refracted = actually incident ray), and writing 
Sj^s for the original incident ray, the finally emergent ray will be 
given by s'-s+SgA, (IV) 

the sum to be extended over all surfaces. 

In the case when all the surfaces are reflectors the last equation 
assumes a remarkably simple form, as will be seen in a later section, 
and offers as a final formula conspicuous advantages. In the 
case, however, of refracting or mixed systems (refracto-reflectors) 
we shall not actually require the resultant (IV), but shall employ 
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the formula (II) successively at each surface, after having developed 
an equally convenient formula for the " transfer " of the data 
concerning the last refracted ray to the next surface. This latter, 
the transfer formula (which will replace and supersede * the formid- 
able array of what the German opticians call " die Uebergangs- 
formeln "), is an unavoidable supplement of (II), the refraction 
formula. The former will be soon developed in one of the fol- 
lowing sections, again without the " help " (which is rather a 
nuisance or obstacle) of any scaffolding or framework of reference. 

2. The Nature of the Data defining an incident Bay. 

It will be recalled from the rudiments of vector algebra that all 
of our vectors, such as r or s, are free or non-localized vectors. 
That is to say, any line segment of length 5 and concurrently 
parallel to a vector s drawn concretely somewhere, is, by the very 
definition of the equality of two vectors, for all purposes of our 
vector algebra, equal to s. In other words, to give s means : to give 
p. and the direction of the ray only. Now, if the refracting surface 
is a plane, this is all that is needed to find the refracted ray s' 
(//, of course, being given). For the whole planet has but one 
normal n. Thus, unless we desire to know, for instance, the glass 
thickness traversed by a ray in a prism, the localization of the 
incident ray s is immaterial and not needed. Not so, however, when 
the refracting interface of the two media is a curved, say a spherical, 
surface. For then n is a function of position on the surface. Thus, 
s giving the direction of the incident ray, we require further data 
for its localization. Now such a datum is obviously the place 
on the surface where it is struck by the ray s. This amounts to two 
scalar data defining position on the surface, let us say u, v, measured 
along any network of curvilinear coordinate axes spread out over 
the surface, as for instance, its lines of curvature. We need by 
no means specify these u, v. Clearly, " to give the incident ray " 
means to give s J and the place (2*, v) of its incidence upon the 

* In the sense that it will be seen to be applicable to cases to which the usual 
formulae are not. 

t Euclidean space being, of course, tacitly assumed throughout. 

: Any two scalar data, the setter of the problem prefers. \Ve do not criunt n 
involved in a=/*r. In general we shall not count data describing the optical 
system, 7't's. refractive indices, curvatures, and so on. 
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surface in question. But this being given, the normal n is also 
given, which being a unit vector amounts precisely to two 
scalar data. 

Thus, whatever the nature of the numerical data for the incident 
ray, be they angles or projections upon any lines or planes, or 
distances from axes (artificial or intrinsic), we are fully entitled 
to say that " incident ray given " means, for a surface <r K) 
S K and n K given. 

Such being the case, formula (II) finds the refracted ray s/. 
Moreover, it finds it together with its localization ; for the latter 
ray starts from the point of the surface o\ struck by the incident ray. 
But in order to apply (II) to the next surface 0-*+ 1 we must transform 
the last acquired data s/ =8*4-1, n 
into 8*4-1, 11*4.1, 

that is to say, we must learn how to replace n* by n*+i . It is this 
which will be given by the transfer formula. Manifestly, having 
once obtained 8*4.1, 11*4.1, we stand before the surface OVH precisely 
in the same position and with the same equipment as before the 
preceding surface. Thus the process of tracing the ray can be 
continued indefinitely. 

There is no essential difficulty in developing, in correspondingly 
general terms, the transfer formula for a pair of surfaces of any 
form. But for all purposes actually met with in technical optics, 
it will be enough to develop the formula for spherical surfaces (of 
any finite radii), and later on (in a special section on prismatic 
systems) also for the sub-case of two planes or a plane and a sphere. 
In the more important case of spheres proper the transfer formula 
is indispensable for finding even only the direction of the refracted 
ray, s'^.i, while in the case of planes (prisms) the ray can be traced 
through the whole system, as far as its direction is concerned, 
without the " transfer formula " ; the latter being required only 
when one is interested in the thickness of the several media 
traversed. 

3. The Transfer Formula for Spherical Surfaces. 

It was said in what precedes that, the point (u, v) on a surface 
being given, the incidence normal n is thereby given. In general, 
however, the converse manifestly does not hold ; two or more or, 
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as in the case of a plane surface, even oo 2 points having parallel 
normals. But if the refracting surface is a sphere, of centre 0, then 
we have a one-to-one correspondence of points (u, v) and normals n, 
there being no two points to which the same n belongs. This 
offers the advantageous possibility of considering the incident 
ray to be given directly by s and n. It will be convenient in the 
case of any spherical surface to draw all the normals n as unit 
vectors from the sphere's centre as their common origin. The 
two vectors s, n being given, we know already, by (II), how to 
trace the ray through the corresponding surface, i.e. how to find 
the refracted ray s'. 

What we still require is to make this ray ready for tracing it 
through the next surface of the system. This is an easy geometrical 
problem, although the optical authorities have made it a compli- 
cated one by treating it with too many " auxiliaries." Treating 
it (according to our principle) intrinsically, it is, and will remain, 
a simple problem. 

Let <TJ, <r 2 , as representatives of any pair cr K , o^+i, be two con- 
secutive refracting spherical surfaces of the given optical system. 
Let O l9 0% be their centres, and R lt R 2 their radii. Draw all the 
unit normals ^from O lt and all n^from 2 , as a rule. The optical 
system, as always, being given, so are lf 2l and all the following 
centres, 3 , etc., collinear or not (centred or not centred system). 
Thus let D be the scalar distance of the two centres contemplated 
and a an unit vector drawn from 0^ towards 2 , as in Fig. 3. (The 




FIG. 3. 

reader should note that the surfaces need not turn their convex 
sides towards one another, nor need X precede 2 ; the formulae 
to be given presently will hold always.) Let HJ be the particular 
normal pointing to P lf the place struck by the incident ray s t 
which penetrates into the second medium as &i' =B 2 , and strikes the 
next surface in P 2 a point to which the particular normal n t 
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corresponds. Having found s/^Sj by (II), we wish to find P 2 , 
i.e. the normal n 2 . 

In order to obtain it walk round the quadrilateral 0i/VV?2 
back to the starting point O v This being a closed polygon, no 
matter whether plane or not, or whether you recross your path * 
or not, the vector sum of the four vectors traversed is nil, i.e., 
calling L the (scalar) length PiP 2 , or medium " thickness " traversed 
by the ray, and remembering that S 1 / =s 2 =/i 2 r 2 , and therefore, 



(i) 

This, in fact, is already the required transfer formula, expressing 
the unknown n 2 in terms of things given all but the scalar co- 
efficient L, a length. In order to find this (which, by the way, is 
not only needed for " tracing " the ray but will serve the computer 
in his further task of calculating " astigmatism and curvature "), 
eliminate from (i) the principal unknown n 2 by squaring both sides 
of this vector equation. In doing so remember that n 2 2 = l, as 
are also the squares of the remaining three Clarendons. Then the 
result will be 



a quadratic equation for L, with 

P-Da-ie^-D-B^ (3) 

which is a given vector. This finds L and completes the required 
transfer formula (i) which can be written, more shortly, 



Is not this simple ? In order to estimate the simplification of 
the deduction, due to the vector method, read the same subject, 
for instance, in the well-known collective work of the scientific 
Jena staff, t where it covers ten pages of print (in which, to be 
just, not one but three alternative " Uebergangs "-methods are 
described, none inferior to the others either in rigour or in clumsi- 
ness). In harmony with this, their " refraction-formulae " are 
equally complicated, the whole apparatus consisting of something 
like fifteen or sixteen formulae. The complication in the scalar 

* Which would happen with the choice of another exemplifying figure. 

t Die Bilderzeugung in optischcn Instrument^ edited by M. von Rohr, Berlin, 
J. Springer, 1904, pp. 52-62* It is, nevertheless, a valuable work. 
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treatment is of course due to the circumstance that skew rays 
are aimed at, not cutting 0^0^ But as the reader will have 
noticed, we have nowhere assumed that the incident ray has any 
such privileged orientation. To repeat, the quadrilateral Oi/VV^ 
need not be plane ; the resulting formulae (i) or (i'), and (2) 
with (3) are perfectly general. 

The solution of the quadratic equation (2) will, in numerical 
cases, cost but little trouble, and this will be repaid by obtaining 
L not only as an auxiliary for the transfer formula, but also as 
one of the " elements " needed for estimating astigmatic and 
associated defects of the optical system in question. The equation 
(2) has, of course, in general two roots, 

L = ( 



corresponding to the two points, P 2 and P 2 ' (f Fig. 3), in which 
the ray cuts the second sphere, and in practice there will be no 
difficulty in discriminating between the root or point (P 2 in the 
case of Fig. 3) to be retained and that to be rejected. And if 
both roots coincide or are complex (conjugate), there is an end 
to the ray-tracing ; for in these cases the ray does not penetrate 
into the next medium. 

It would be entirely superfluous to enter at this stage into the 
details depending on the configuration of O lt O f , and the following 
centres, together with the sense of the curvatures (concave, convex ; 
convex, convex, and so on). The reader will do well to draw several 
figures for the possible cases, especially for those in which the sides 
of the quadrilateral cross one another, and to inspect them care- 
fully. Among other things he will then convince himself that, with 
the vector a drawn always from Oi to t (whether O l precedes O, or 
follows it), formula (r), supplemented by (2), (3), remains always 
valid. 

Having obtained n 2 by means of the above formulae, w$ can 
proceed to trace the ray through the surface having its centre 
at 2 ; similarly we shall pass from that to the surface centred 
at 8 , and so forth. No complications will be met with if it should 
happen that the system is not centred, i.e. if O lt 2t 0$, 4 , etc., 
are not collinear, or even not coplanar. For the first two surfaces, 
a=a 12 will be the " axial " unit vector ; for the second and third 
another unit vector &%$, drawn from 2 towards 3 , will replace 
it, and so on. These vectors, being not artificial but intrinsic, 
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inherent in the optical system, cannot give rise to any superfluous 
complications of the calculator's work. If the centres happen 
to be arrayed in a straight line, all the better : in that case a single 
unit vector a will perform the function of the axis throughout 
the tracing process. The calculations will become correspondingly 
simpler, again for an intrinsic reason, to wit, because the optical 
system as such is a particularly simple sub-case of the general 
non-centred system. 

4. Summary of the Tracing Process, for Systems of 
Spherical Surfaces. 

It will be well to summarize the whole procedure of tracing a ray 
through any such system. 

The incident ray being given in direction and position, the data 
of the problem are, in whatever form, 

B! and n r 

From these, /i x cos t^-s^ and t/ by /*/ sin h'-^i sm *'i follow 
immediately, whence the scalar coefficient 

ft^x'cosiy-MiCOSiV* 
This being found we have, by the refraction formula (II), 



which is also B 2 , the ray incident upon the second surface of the 
system. 

The vector drawn from O l to the next centre 2 being D 12 
and 



the medium thickness L is calculated as the appropriate root 
L - (Pr 2 ) 



of the quadratic equation (2), with r a 

Finally, by the transfer formula (i) or (i'), 

n, -|J- (&*-*) 

^2 

Thus the two vectors are found, 

s 2 and n 2 , 
ready to trace the ray, thus localized, through the second and 

* Which is to be replaced by - 2/*i cos i lf when the surface is contemplated as a 
reflector. 
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to carry it over to the third surface, and so forth, up to the last 
surface of the system, giving the finally emergent ray through 
two vectors, say, 

*' and n', 

the first determining* its direction, and the second localizing it 
completely with respect to the system. 

The results, at any intermediate surface o- K , as well as at the 
last surface of the system, are all of the same form as the original 
data (B V n x ). 

Each of these pairs, whether of data or results, amounts to the 
knowledge of four scalar numbers, two involved in s (not counting 
p.) and two in the unit vector n. 

This exhausts the matter in its general aspect, the procedure 
being valid for any orientation and position of the incident ray 
and for any system of spherical surfaces, centred or not. 

With regard to its particular or practical aspect, the successful 
application of this compact set of vector formulae to numerical 
cases will depend only upon one negative condition, which is : not 
to spoil their intrinsic simplicity by inappropriate splitting or 
decomposition of the vectors. A large portion of the remainder 
of the present work will be dedicated to an attempt to show the 
reader how to avoid this danger. 

As a matter of fact, readers skilled in the handling of vectors 
will, after what has been said, hardly require any further explana- 
tions. This little book, however, being chiefly intended for those 
who are not experts in the use of vectors, it has seemed indispensable 
to add some hints at least which are likely to help the majority 
of readers in applying the general formulae to more or less concrete 
cases. These will be given in some of the following sections. 

5. Hints and Examples. 

Having hitherto, in developing the fundamental and general 
part of the subject, followed a systematic course, we can henceforth 
without danger of confusion depart therefrom and choose our 
subjects freely. Our aim being to explain the general procedure 
by illustrations and hints, a pedantic observance of order would 
only hamper us unnecessarily. 

Let us start with the simple case of a lens consisting of a single 
piece of glass of index /A limited by two spherical surfaces or,, <r s , 
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the surrounding medium being air or empty space. To fix our 
ideas let us take a 

Biconcave Lens. Let, as before, O lt RI be the centre and the 
radius of the first surface, upon which the incident ray s = r 
impinges, and 2 ^2 t ^ e centre and the radius of o- a , the back 
surface ; finally let B lf n 2 be the unit normals, drawn from the 
centres towards the points at which the ray enters and emerges, 
as in Fig. 4, which however is scarcely needed. By assumption, 




FIG. 4 . 

the two unit vectors r, n are given, amounting each to two 
numerical data. The first practical question is : How are they 
given ? Generally speaking, this will depend upon the decision 
of the person who sets the problem. Yet, since we are here setting 
our own problems, let us put the data in a reasonable form, seeking 
only to avoid artificial frameworks of reference. Now, the lens, 
consisting but of two spherical surfaces, offers us but one intrinsic 
directed element, viz. its axis, which we will represent by the unit 
vector a drawn from O l towards 2 , these points themselves 
having already been exploited by taking them as the origins of 
the normals. Besides a there is nothing intrinsic to build with 
a framework.* Now, if the ray r happens to lie in a meridian 
plane (passing through the axis), then nothing more than a is 
needed for the ray's relevant definition ; say the angles r and n t 
make with the axis. But let r be a skew ray, not contained in any 
meridian plane.t Then the lens itself offers us manifestly nothing 
to complete the determination of the direction r. But the ray 

* A plane normal to a may be a welcome auxiliary, but gives, manifestly, nothing 
independent of a itself. Nor shall we look for the rim of the lens to help us. 

t The simpler case of a meridian or non-skew ray may be taken up as a sub-case, 
and as an easy but useful exercise. 
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itself, as localized, does offer us something, and exactly as much 
as is needed. In fact, it strikes the first surface in a given point P l9 
given in one way or another. Through that point draw your 
41 first meridian " plane, marking it with any unit vector b, say 
perpendicular to a, and count, from that plane or semiplane, the 
geographical longitude <, the companion of the polar distance 
or colatitude tf^^a, n x . Thus 0, <f> will be a polar system of 
coordinates with centre O v Remember that we have drawn our 
first meridian plane a, lf b, only through the end point P l of r, 
but that the ray is in general not contained in that plane. 

Now, this is all we want. The vector a being furnished by the 
lens, and b by the incident ray, we have everything which is 
required to write down numerically (scalarly) the data. In fact, 
" r given " means now ra, rb given, say : 

ra = cosa, rb = cos/3. (5) 

For symmetry, although it is superfluous, introduce a third unit 
vector c = Vab (normal to a, b, towards the reader), and write 

re = cos y, so that (ra) 2 + (rb) 2 + (re) 2 = I . (6) 

Keep in mind, however, that the first two are by themselves suffi- 
cient. This settles the question of r. Next, "nj given" now 
means lt ^ given, i.e., to be general, 

n x = a cos l + (b cos fa + c sin <#> x ) sin O l9 (7) 

and since, by the choice of the first meridian, <i=0, 

n x = a cos O l -f b sin lt (70) 

In fine, the complete numerical data, immediately attackable 
by the computer, are 

lt <h-o, ra, rb, 

the latter two as in (5), or, if he prefers, with re as in (6) instead 
of rb. 

Such being the scalar data, act according to the compact pro* 
cedure of Section 4. Thus * 

cos ij = r^ = cos l (r x a) -f sin t (r x b) 

= cos #! , cos a + sin l cos /3, by (70), (5) ; 

whence sin i^ - sin i t and g l = p. cos i/ - cos ^ , 

* Remembering that in the present case Bj=:ri=r, ^= 
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ready for the first refracted ray 



If the reader wishes to split this at the present stage he can do so 
at once, multiplying (8) scalarly by a, b or c. But, being interested 
only in the finally emergent ray, let us proceed according to the 
general prescription of Section 4. To " transfer " the ray r a , 
found by (8), to the surface <r 2 , write D for the distance O l 2t and 
P = .Da - Rin lt as before ; then the scalar product required for L 
will be 



), by (8), (5), (?a), 

and, the square of P, 

/* = /> + /?!*- 2DR l cos 9 l9 by (7 a). 

This finds the glass thickness traversed, which in the present 
case is, manifestly, the smaller of the two roots of (2), 

L = (Pr 2 ) - v/(S^P* + ., 

which the reader will properly contract by using the last two 
expressions. The " transfer " is thus completed, giving n 2 , and 
therefore also the point P 2 , where the ray strikes the back surface, 



To know the numerical coordinates of this point may be of 
immediate interest. Now, calling 2 , < 2 its colatitude and its 
longitude, with 2 as centre, but still using the previous first 
meridian plane, in short the same unit vector b, we have just, 
as in (7), 

n 2 = - a cos 2 + (b cos <f> 2 + c sin < 2 ) sin 2 , 

i.e. for the required two angular coordinates, 

COS 2 = ~ n 2 a COS ^2 S * n ^2 ~ n 2^i S ^ n ^2 S * n ^2 = n 2 C > 

or, ultimately, cn 

cos 2 - - an 2 ; tan < 2 = ^- 2 . (9) 

Notice that ^ being zero, </> 2 need not necessarily be so; in 
other words, n 2 is in general not coplanar with a, n r It may be 
well to recall here that, for any three vectors, 
AVBO=BVOAOVAB 
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is the volume of the parallelepipedon constructed on A, B, as 
edges, and that, therefore, the sufficient and necessary condition 
of coplanarity of n lf a, n 2 is the vanishing of the (scalar) expression 



or, since P is itself coplanar with a, n lf of the expression 

/iaVn 1 r 2 =aVn 1 r 1 , by (8). 

We thus see that < 2 ^o, ** that P 2 is outside the meridian plane, 
unless the incident ray r x happens to be in a meridian plane. For 
any skew incident ray we have < 2 = o. 

The coordinates </> 2 , 2 of the leaving point P 2 being thus found 
by (9), the emergent ray is localized, and it remains only to find 
its direction r 2 '. This is, by the refraction formula, 

*i =f*r 1 +g t n tl 2 =cosiy-/*cos/ 2 , (10) 

the latter to be found at once from the known 
r2n 2 =~(L-Pr 2 ), 

as was g, a moment ago. This completes the tracing of the ray. 
For, if a', /?' be the angles which the finally emergent ray makes 
with a, b, we have only to multiply (10) scalarly by a, b, to obtain 

cos a! = r a 'a, cos /?' = r 2 'b. 
Thus the original incidence data having been put in terms of 

0i, </>i=o, a, ft, 
the final answer is obtained in similar language, viz. 

<**, * ', p. 

This essentially exhausts the problem of a simple lens. If 
instead of a biconcave we have a biconvex, or a " vex-cave " lens, 
etc., the only modifications needed will concern the sign of a 
in the expansion of n 2 , and the sign of the square root appearing 
in the " glass thickness " L traversed. The exhaustive discrimina- 
tion of all possible cases, which will among other things depend also 
upon the relation of R lt R 2 , D, may be left as an easy, but useful, 
exercise for the reader. The central glass thickness is, in the 
above case, D"(R l + R 2 ) f and similarly, mutatis mutandis, for the 
remaining cases. 

Any further information concerning the emergent ray may, 
whenever required, be easily deduced from the obtained 2 , </> 2 ; <*', /?', 
or vectorially, from n 2 , r 2 '. Thus, for instance, to decide whether 
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the emergent ray meets at all the axis ot the lens (at a finite or 
infinite distance), test the triple product 

/>=-aVr 8 'n 2 , (n) 

which is a scalar, or ordinary number. If p=o, the emergent 
ray or its prolongation meets the axis somewhere or is parallel 
to it ; and if f = 0, then r 2 ' is a skew ray. By (10) and (8) we have 

p=aVn 2 (r 1 +^ 1 n 1 ), and so on. 

Let the reader complete the development of p, and discuss the 
circumstances of its vanishing or not. A useful exercise. 

As a further exercise, assume (abstractly or in a numerical example) 
that/> = o, so that the emergent ray is a skew ray, and find, vectorially, 
by using r t ', n t , the shortest distance of the emergent ray from the 
lens axis ; also the point of the axis which lies at that distance 
from the ray. 

As another exercise treat the trivial case of p=o and find, again 
vectorially, the distance, from the back surface, of the point of 
intersection of the ray r t ' and the axis. 

Construct a handy formula for the angle between the original 
incident ray and the finally emergent one, remembering that its 
cosine is lyr/. 

The reader will also find it useful to invent for himself similar 
problems, resolve them vectorially and evaluate the appropriately 
scalarized expressions numerically. In this way only will he acquire 
the needed skill in handling not only vectoriai optics but vector 
investigations in general, which, after all, will be found also to be a 
comparatively pleasant occupation. 

Not confining himself to lenses, let him also treat a single or two 
spherical mirrors, for which g^ assumes the simpler form - 2 cos i K , 
as explained before. 

In actual numerical ray tracings the principal and the accessorial 
vector formulae will best be written out (and glanced at whenever 
necessary) in the margin of the computer's sheet, or on an independent 
little sheet of paper conveniently placed before his eyes. Now and 
then it will be found convenient to note a vector result by writing 
out the figures found for its " components/ 1 as for example 



having previously decided upon the choice of the length unit. It 
will be remembered, of course, that every unit vector is dimensionless, 
being a vector divided by its own (scalar) size or representative 
length. Finally, details concerning the manner of writing and 
arranging the numbers or their logarithms will best be left to the 
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individuality of each reader. It is a familiar experience that, within 
very wide limits, everybody finds his best way of arranging the 
numbers by himself, after he has repeatedly executed a given type 
of computation. It is by no means desirable to drill away human 
individuality. 

6. Prismatic Systems. 

Leaving for the present the spherical separation surfaces, let us 
turn to plane ones. A series of different optical media separated 
from one another by plane surfaces a-^ constitutes a prism or 
prismatic system. This is of almost equal interest and importance 
with lens systems. As a matter of fact, any prism system is but 
the simplest sub-case of what has been treated in the first two 
sections. Yet, in order not to omit any opportunity of familiarizing 
the reader with the vector treatment of optical problems, let us enter 
somewhat into the details of these comparatively simple systems. 

A single Prism. Let us begin with the case of but two plane 
surfaces er lf <r 2 limiting (to an extent) a medium, say glass, of 
refractive index p.. Outside /A = I . In the first place, let us treat 
the rays apart from their localization, i.e. only as far as their 
directions are concerned. Each refracting face is then fully 
characterized by a single unit normal, n t for the first, n 2 for the 
second plane face. Remembering that inversion of the n's leaves 
all the fundamental formulae intact, draw both normals towards 
the air, as in Fig. 5, so that the prism angle or refracting angle 
is given by 

n 1 n 2 = cos(7r-)= -cose. (12) 

The question of an " auxiliary " framework now entirely dis- 
appears. Unlike the simple lens {or any centred lens system) the 
prism is liberal enough to furnish us, by its own intrinsic elements, 
a complete natural reference system. For it has, together with n lf n 2 , 
a well-determined edge. Take a unit vector e along the edge, 
normally to the plane of Fig. 5, and towards the page (or down- 
wards if your reading table is horizontal), so that 

Vn x n 2 =e sin (w - e) =e sin . (13) 

The angle being given, of course, this is a simple relation between 
the three natural reference vectors. Returning to Fig. 5, drop 
the outline of the prism itself, shift the normals n lt n t to a common 
origin, and draw instead the unit edge e, as in Fig. 5a. The prism 
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will then be represented by the simple three-vector skeleton, with 
fji subintended. 

Now let any ray r impinge upon (the " first " face of) the prism 
thus represented. This given ray may also be drawn from the 





FIG. 5. 

same origin, though it is scarcely necessary. It finally emerges 
from the prism as r a '=r', say, viz. by (IV), Section I, 

*'=*+gi*i+gP* (14) 

Here r is, of course, in general any skew ray. The explicit 
expressions for the scalar coefficients g lt g 2 can be found at once. 
But even without the knowledge of their values the vector formula 
(14) shows us at a glance an interesting and perfectly general 
property of the prism. Multiply it scalarly by e ; then the result 
will be 

r 'e = re = cos 0, say. ( 1 5) 

Thus the final emergent ray makes with the edge the same angle 
as the original incident ray, whatever its orientation. In other 
words, the only effect of the prism is to turn the ray rigidly round 
the prism edge as axis. The amount of this rotation depends upon 
the properties of the prism, n and c, and upon the orientation of 
the incident ray. 

The ray r being given, so is the incidence angle i=sarccos (rax), 
whence also the coefficient g l9 appearing in the resultant formula 

(14), 

gj = v/x 2 - sin 2 i - cos i. ( 1 6) 

Similarly the second coefficient g 2 could be written down, as in the 
case of the lens (and somewhat easier). But it is more elegant to 
represent fully the final emergent ray r' without the help of # 2 .* 

* Especially as the orderly way of calculating g would presuppose the know- 
ledge of the final exit angle t'=/ 2 '. 
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To get rid of this troublesome coefficient multiply (14) vectorially 
by n 2 ; then, by (13), and remembering that Vn^^O identically, 



. 

This expresses a vectorial property of the finally emergent ray 
in terms of things which are all known. Together with the capital 
property (15) the last equation completely determines the emergent 
ray in its orientation. 

Exercises. Discuss equation (17), deducing from it various 
vectorial consequences and translating each of these into geometrical 
language. 

Find the expression for the amount of rotation round prism edge 
as axis, which converts the incident r into the emergent r'. 

Find, in the general case of any skew r, the limit condition, beyond 
which the ray does not emerge but is internally reflected (at or t ). 

Treat vectorially the trivial, but instructive, case of r _L e. (Assume 
at the outset r of the form an t +6n t ; 0, b scalars.) Translate your 
results into scalar language, deriving thus the usual formulae. 

In practice (spectroscopy, etc.) the aim is, of course, to make 
re=o, or rj_e. Yet the manufacturer, and the student in his 
physical laboratory, can but approach this ideal within, say, 10" 
or even broader limits. It is, therefore, important to estimate 
the effect of this error upon the emergent r'. Develop, therefore, 
the above formulae up to the very end by assuming a value of 



slightly differing from zero, i.e. equal to a small fraction. In other 
words, put 0=-;t8, where 8 is a small angle, say 10* or i' or 5'. 
Employ, whenever necessary, the well-known developments in series : 




Find the effect upon r' of a small change or " variation " in the 
direction of the incident ray r. Replace r by r+8r, remembering 
that r l =i, and therefore, rSr=o, or SrJLr, and so on. 

Similarly study the effect of a slight variation Sc in the refracting 
angle of the prism, as expressed by (12). Scalarize the resulting 
formulae and interpret them in physical (geometrical) language. 

Medium Thickness traversed. In what precedes, the rays 
incident upon, transient through and emerging from the prism 
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were considered with regard to their direction only, independently 
of their localization. In fact, it Was possible 'for that purpose to 
replace the prism by a mere triad of unit vectors, n lf n a , e. The 
" transfer formula " was not needed owing to the circumstance 
that the same n a served as incidence normal for the whole second 
plane. This is the reason why the treatment of a prism is con- 
spicuously simpler than that of a lens. In short, the tracing of 
rays through a prism as far as their direction goes requires no 
" transfer " process. Some device of that kind, however, becomes 
necessary if we desire to find the length of the ray, whose direction 
is ^'=12, within the prism, or, as we will call it for brevity, the 
medium thickness traversed. Let this be L. 




FIG. 6. 

Obviously, in order to find it, we must localize the incident ray r 
with respect to the prism, by giving, say, the point P l where it 
strikes the first face <r 1 . Since P cannot be defined by % (as in 
the case of a spherical interface), the normal being common to all 
points of the plane, it is indispensable to determine position on cr 1 
by some coordinates. As such we may conveniently take y lf 
measured along e (the edge), from an arbitrary origin on the 
edge, and x l perpendicularly to the edge. Similarly y a > #2 f r ^ e 
second face cr 2 of the prism. We can collect these pairs of co- 
ordinates by introducing the vectors 

(19) 



where kj, k, are two unit vectors along the axes of x ly # t , as in 
Fig. 6. 

To obtain L, walk round the triangle OP^O ; then the result 
will be 

(20) 
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This, in the present case, is the " transfer formula," somewhat 
simpler than in the case of spherical surfaces. Multiply (20) scalarly 
by n 2 , which is normal to p 2 ; then 

Lrjjn^-p^, (21) 

which finds L, and completes at the same time the prismatic 
transfer formula. Also, of course, L 2 = (p t -P 2 ) 2 . 

Exercises. In (21) substitute r a =ri', using the fundamental 
refraction formula, thus completing the little investigation. 

Study L in its dependence upon the prism angle , and upon the 
direction and the position of the original incident ray. 

Deduce the formulae for the simplest sub-case re=o. 

Find the deviation of the ray, noting that its cosine is rr'. Establish 
the conditions of minimum deviation. 

System of any Number of Prisms. The tracing of any ray, 
in direction and position, by the successive application of the 
refraction formula _ _ / 



and of the transfer formula typified by (20) with (21), offers no 
difficulties, nor essentially novel points, and need not, therefore, 
detain us here. We can imagine that by the application of the said 
process all the scalar coefficients g* are already determined by the 
given incident ray r and by the properties of the given system 
of prisms. 

Thus, supposing, for instance, the first and the last medium to 
be air, the emergent ray r' will be, by (IV), 

r' =r + 2&II* =r +g l n 1 + ... +g lt n lc , (22) 

where n* is the unit normal of Kth plane face. In this more 
general case (* >2) it will be convenient to draw all the normals 
n lt n 2 , n s and so on in one and the same sense, say, from the first 
towards the second, from the second towards the third medium, 
and so on. Again, as in the case of a single prism) the whole 
prismatic system can be represented by a skeleton of coinitial 
unit normals n lf n 2 , n 3 , etc., and unit edges e 12 , e^, etc. These 
latter, no matter whether coinciding in space or not, can also be 
shifted to the same origin. (An air gap will also count as a com- 
ponent "prism," with fii.) Thus a system of AT-i prisms, 
i.e. of N interfaces, will be entirely represented by a skeleton of 
N vectors n K and N - 1 edge vectors e 12 , etc. If e lf c 2 , etc., are the 
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prism angles, the relationship of the latter to the former unit 
vectors will be expressed by 

sin c t . e ia = Vn^j, sin < 2 e 23 = Vn^, etc. 

Any given incident ray r impinging upon this system will be 
converted into r', as expressed by (22), where the g K are assumed 
to have already been found. (The system itself constitutes, of 
course, the best and more than complete framework of reference, 
thus giving no place for artificial " auxiliaries.") All properties of 
the emergent ray can at once be derived from the vector equation 
(22). We know already that the sole effect of each prism taken 
separately is to turn the incident ray round the edge as axis, 
leaving the angle with the edge intact. The effect of the whole 
system will be the resultant of a number of such rotations round 
e 12 , 623, etc., as axes, taken in their proper order, of course, if, as 
in general will be the case, the edge vectors do not coincide, i.e. if 
the prism edges are not all parallel to one another. 

The optical manufacturer attempts in ordinary circumstances 
(as in the case of spectroscopes) to make the prism edges parallel, 
i.e. e 12 =e23=e34, etc. But, of course, he can succeed in doing so 
only approximately. Let us call such an ideal system a correct 
one, and any other, with edges more or less deviating from 
parallelism, a defective prism system.* Now, especially if a system 
of prisms is defective, formula (22) may be very helpful, because 
with n lf n 2 , n 3 , etc., not coplanar, any of the ordinary Cartesian 
methods would soon lead to the greatest confusion and would 
require very careful drawings and possibly three-dimensional 
models to follow the rays in their tortuous course. The vector 
formula will therefore be particularly useful if it is required to 
estimate the effect of, say, small deviations from parallelism upon 
the emergent ray. 

If the system of prisms happens to be (or is considered as) a 
correct one r then all the normals n* are coplanar and all the unit 
edge vectors are fused into one, e, say. Multiply (22) by this 
vector ; then the result will be 

r'e=re, 

i.e. 0' = 0, precisely as for a single prism, as could be expected. 
Whatever the orientation of the incident ray r, the prism-system 

* Although the reader may well imagine cases in which non-parallelism is a 
feature deliberately aimed at. 
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turns it rigidly round the common representative edge as axis, 
leaving its inclination intact. Notice also, that all the normals 
n,c being coplanar, the vector sum in (22) is again a vector in their 
plane. Such a vector can always be represented as the sum of 
two vectors along two fixed unit vectors in (or parallel to) that 
plane. Let these be a, b. Then (22) will assume the form 

r'=r + ^a + b, (23) 

where A, B are some scalar functions of the orientation of r, and 
of the characteristics of the prism system, i.e. of the p, K and the 
relative orientation of the n K (the latter already including all the 
angles C K ). These scalar functions A, B can easily be determined 
for the case of any concretely prescribed prism system. The 
resultant deviation (angle) ^ = arccos (rr') follows at once by 
multiplying (23) scalarly by r, which gives 

cos ^ = i + A cos (r, a) + B cos (r, b). (24) 

As a useful exercise the reader may treat fully the case of two 
glass prisms separated by an air-gap, either as an " ideal ".or as 
a " defective " system. 

In the latter case it will be of particular interest to fully develop 
the formula for the effect of a small deviation from edge parallelism. 

7. System of Plane Mirrors. 

Let again, as in Section 6, each interface <r K be plane, but let 
each of them be contemplated as a reflector. If the position of 
rays is disregarded, such a system, which may shortly be called 
a multiple reflector, will be entirely determined by the assemblage, 
and the order of succession, of the unit normals n lt n 2 , n 3 , etc., 
drawn, let us say, away from the reflecting faces. 

As was already hinted at on p. 3, the resultant formula (IV) 
assumes for any such multiple reflector a particularly simple 
form. 

In fact, returning to the reflection formula (I), written down 
at the very beginning, and using now a single dash for every 
reflected ray, we have, for the first reflected ray, 

r 1 / r 1 -2n 1 (n 1 r 1 ), (25) 

where r x is the (original) incident ray. We already know that 
this is a particular case of the refraction formula (II), page 2, 
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viz. with /x' replaced by -/*. But there is this important difference, 
that while (II) contained in its right-hand member, through g, 
the " unknown " s', formula (25) contains in its right-hand member 
T I only (besides the given n x ) and no trace of r/. It is true that, 
as was promised on p. 2, the vector s' can easily be eliminated 
from the right-hand member ; but in doing so we obtain for the 
refracted ray a formula which in its structure differs essentially 
from the reflection formula (25). In fact, squaring (II), and 
remembering that s /2 =fi' 2 , s 2 =/x 2 , we have 
g 2 + 2(sn)g=/* 2 -/i 2 , 
whence g, in terms of sn =//- cos z, and by substitution in (II), 

r ' = P T - /o!n (nr) + nV^'+'p^far) 2 , (26) 

where /o= /*//*' and <r 2 = i-p 2 . Here r only appears in the right- 
hand member ; but the operation to be performed upon r in order 
to obtain the refracted ray r' is * manifestly a wan-distributive one. 
That is to say, if (26) is written for the moment 

r' = o>r, 

w standing for an operator, then w(A+B) is not equal o>A + wB. 
It is for this reason that we have said (Section i) that it is not 
worth while to eliminate s' from g. Nor would it give us any 
advantage to study the complicated and non-distributive operator 
cu, as defined by (26), which might be called the " refractor." 

On the other hand, a glance at the reflection formula (25) will 
suffice to show that the operator converting the incident into the 
reflected ray is distributive. This is the reason of the simplicity of 
multiple (plane) reflection, corresponding to a repeated application 
of such operators. 

In fact, nj(A+B) is xijA+i^B, which proves the distributivity 
of the operator fl lf say, involved in (25). That is to say, if we 
write that equation r x ' ^ = 12^, we have, for any A, B, 



This is a capital property of the " reflector " 12, as contrasted with 
the non-distributive property of the '* refractor 1 ' w. 

In order to see the structure of the 4< reflector/* belonging to the 
mirror whose normal is n ly write (25), using the dot (instead of 
brackets) as separator, 

r 1 '=r 1 -2n 1 .n 1 r 1 . 

* Quite apart from being more complicate than the reflecting operation. 

c 
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This is but a slight change in notation. Next, instead of writing 
twice the same operand r lf put 

r 1 '=[i-2n 1 .njr 1 , 

to mean, by convention, precisely the same thing as before. Thus, 
the operator Q tt converting any incident ray into the corresponding 
reflected ray, is seen to be 

12 1 = i-2n 1 .n 1 , . (27) 

n t being the mirror's unit normal. Such an operator is called a 
dyadic (n l .n 1 being called a dyad, since it contains two vectors 
separated by a dot). The reader need not, however, be deterred 
by new names or new concepts.* Since he knows how fij was 
introduced, there can be no mystery about the meaning or the 
properties of this operator. In fact, he can read the definition (27) 
in plain English thus : to operate with 12 X upon any given vector A 
means to cut off twice its projection upon n t along n t and to subtract 
the vector thus obtained from I . A, that is, from A itself. This 
is a very simple operation indeed, more simple than most of those 
the reader daily performs, whether intellectually or physically, 
There will thus be no danger of a misunderstanding in using this 
operator, which is but a very particular example of the powerful 
dyadic. 
Thus our first reflected ray will be 

*i'-fli*i. (28) 

Let now the reflected ray r x ' impinge upon the second mirror of the 
system. Then, $2 a = i -2n 2 .n 2 being the "reflector 11 of this 
mirror, the second reflected ray will be, exactly as before, 



and substituting i/ from (28), 

r^-QA'i, (28*) 

which means : operate with i\ upon r t , then with 12 2 upon the 
result of the first operation. This is in essence the same thing as 
to multiply a given number x first by 3 and then the result by 5. 
Both kinds are intelligible and reasonably simple operations. The 
only important difference between them is that while the order 
of the operations in the latter case is a matter of indifference, it 

* I mean, new to the reader. The name and the concepts, due to Gibbs, are 
old enough. 
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is not so in the former case, ft^ is (in general) not the same 
operation as fljftj. So also is the actual reflection from mirror I. 
followed by a reflection from mirror 2. not the same thing as a double 
reflection in the reversed order (unless the two mirrors are in a 
particular relation to one another). 

By (280) it can now be said that the " reflector " of a double 
mirror (n lf n 2 ) is 



Similarly for three, four and any number of plane mirrors. 
The " reflector " of a multiple mirror, consisting of N plane mirrors, 
acting in the order I, 2, ... *, is 



or, more shortly, 

fl =H12 X , where 12* = I - 2n K . n K , (29) 

the " product" extending over all the component mirrors in their 
prescribed order. If any ray r impinges upon the multiple mirror, 
then the final ray, reflected in the said order, is 

r'=12r. (30) 

The problem is thus reduced to building up the resultant operator 
12 from the component ones, 12^. Now, this is done almost as 
simply as the algebraic multiplication of a succession of poly- 
nomials. For to any succession of operators 12, belongs the 
associative property, i.e. 



and similarly for any number of component reflectors.* The 
said simplicity of handling reflectors 12* and of building them 
up into the resultant reflector 12 for the multiple mirror is based 
upon these two properties, the distributivity and the associativity 
of these operators. In each case the process will consist in simply 
*' multiplying " out the dyads n . n contained in the component 
reflectors. The only precaution heeded is to keep the order intact 
and to remember that juxtaposed vectors .not separated by a dot 
are fused into ordinary scalar products. Thus, for example, 

n t . n^g . n 2 =n 1 (n 1 n 2 ) . n 2 =a 12 n 1 . n 2 , 
where a 12 ^BI^ = cos(n lf nj) is an ordinary scalar number. Similarly 



* Henceforth " reflector " will be used for the operator and " mirror " for the 
reflecting surface to which it belongs. 
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for the product of any number of dyads. Such a " product " will 
always be again a dyad, including an ordinary scalar factor. 

In all the following developments we shall employ the short 
notation 

n t n, = cos (n t , n,) =a t ,, (31) 

so that, obviously, 

&IK 

Before passing on to examples, notice that equation (25) or (28), 
when squared, gives 



Thus, the incident r x being a unit vector, so is also the reflected r^. 
The same property holds, of course, for any multiple mirror, 
i.e. if r' = 12r, as in (30), then r' 2 = r 2 = i. Notice also that the 
iteration of any 12, is I or an idem factor ; in fact, 

12,12," i -411* . UK + 4(n, . n,) 2 = I -411, . n, +n, . n, 
(since n, 2 = i), whence 12,12, or 12, 2 = I, which is an obvious property. 
If ^ be the angle contained between the incident and the finally 

reflected ray, we have 

cos^=rflr, (32) 

for any multiple mirror. The expression on the righj hand is the 
scalar product of r into the vector fir, or, which turns out to be 
precisely the same thing, of the vector ri2 into r. 12 being a dyadic, 
composed of dyads a . b, can manifestly be used either as a post- 
factor or a prefactor. For any of the component dyads can be used 
so. Thus 

a . br means a(br) 

and ra . b means b(ar), 

both being perfectly definite results (vectors), in general differing 
from one another, of course. But r(12r) is the same thing as (r!2)r, 
so that in (32) and all similar equations no brackets are needed. 

The reflector 12 belonging to any N-uple mirror will be obtained 
by multiplying out (with the said precautions) all the component 
reflectors 12,. Thus, to* obtain 12 for a quadruple mirror, say, we 
must necessarily have first the operators for a triple, and, before 
that, for a double mirror. Similarly for any N. It is therefore 
but natural to begin with the double mirror, and, before passing 
on to the triple and more complicated mirrors, to stop a while at 
each of the intermediate stages, discussing and interpreting the 
results. It will be understood that by an " Af-uple " not neces- 



SYSTEM OF PLANE MIRRORS 27 

sarily a system of N different simple mirrors is meant, but only that 
N successive reflections (in a prescribed order) are contemplated. 
An incident ray can be reflected from two mirrors, for instance, 
three or very many times, as is the case, for example, in the plan- 
parallel piece of a Lummer plate. 

Double Mirror. The cosine of the angle between the two com- 
ponent mirrors being 12 =n i n 2> as m (3 1 )* we nave , by (29), for 
the reflector of our double mirror, 

12=12 2 fl 1 =[i -anjj.njti-an! . nj, 

i.e., with due regard to the prescriptions of the u multiplication " 
rule, 

ft = i - 2% . B! - 2n 2 . n 2 + 4a 12 n 2 . n lt (33) 

or, introducing the vectojr m =n t - 2<z 12 n 2 , 

12 = i - 2m . nj - 2n 2 . n 2 . (33a) 

This, of course, is valid for any operand, that is to say, for any 
incident ray r, giving r' =12r for the reflected ray. Thus 



where r 1 =rn l , r 2 =rn 2 are the given projections of the incident 
ray upon the mirror normals. Since the edge e, along which the 
two mirrors intersect, is perpendicular to n 1; n 2 , and therefore also 
to m, we have again, as in the case of a prism, 

r'e=re=cos#. 

In fact, the double mirror is but a sub-case of the prism (for /*' 
replaced by -/x). The angle 6 which the ray makes with the edge 
remains intact after the double reflection. This is a general 
property of the double mirror independent of the value of a 12 . 

From (33) we see that 12^ differs from ft 1 12 2 since the last 
dyad n 2 . E^ is not symmetrical (n 2 , n a are two different directions). 
Thus, if a beam of parallel rays r is broad enough to impinge upon 
both mirrors, it is split by the double mirror into two beams, r' as 
before and r" =12 1 12 2 r, so that 

r' -r" ==4^12^2 - i -i - n 2 ]r. (34) 

In particular, if the double mirror is orthogonal, that is to say, 
if a 12 =n 1 n 2 =o, we have r'=r" 

i.e. the two reflected beams coincide. In this case the operator (33) 
becomes a synynetrical (or self -con jugate) dyadic, viz. 
ft =ft 2 ft 1 =ft x ft 2 = i - 2n x . n x - 2n 2 . n 2 . 
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Again, since BI JLn 2l we have 



for if i, j, k, be any normal system of unit vectors (as is n ]} n 2 , e) 
we have manifestly [i .i+j . j +k .k]VV, for any vector V. 
Therefore, for an orthogonal double mirror, whose unit edge is e, 

ft- i+2e.e. (35) 

That is to say, the reflection from such a mirror is, independently 
of the order (12 or 21), wholly equivalent to the reflection from a 
simple mirror whose normal is e, followed by a simple reversal ( - 1). 
This holds for any incident ray r. If the incident ray happens to 
be normal to the edge, i.e. if re=o, we have 

r'=flr = -r, 

that is to say, the ray is sent back parallel to itself. This simple 
property is also familiar from ordinary geometrical constructions. 
But none of the cases of a skew r (oblique to e) will be found so 
easy to treat by the usual methods, based on a good drawing and 
spherical trigonometry, as by the proposed vectorial or dyadical 
method. 

Returning once more to the general formula (33), we have for the 
angle ^ between the incident and the finally reflected ray, by (32), 

i( I - cos vO = rf + r 2 2 - 2a 12 r 1 r t ; 
on the other hand, since r=r 1 n 1 +y 2 n 2 -h(re)e is a unit vector, 

rj 2 + r 2 2 + 2a n r l r 2 = I - (re) 2 = sin 2 0, 
so that the last equation can be written 

|( I - cos \) = sin 2 - 40 12 *V 2 , (36) 

giving the angle ^ = r, r' for any double mirror and for any given 
incident ray r =r 1 n 1 +r,/L 2 . 

From the general formula (36) we see without difficulty that there 
is no such double mirror which would send back (parallel to its 
own path) every incident ray ; in short, that a double mirror cannot 
be a " central n mirror, as are those used in modern signalling. 
In fact, cos ^= - i would*imply 40i 2 *y 2 + (re) 2 =:O, anc * ^ s cannot 
be satisfied for all directions of r. 

Further discussions of the double mirror will be left to the reader 
as an exercise. In this case the adaptation of the vector formulae 
to numerical purposes is almost immediate, and calls? therefore, for 
no further explanations. 
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Triple Mirror. The resultant reflector in this case is Q 
the product of (2 3 = i -2n 3 . n 3 into the operator of the preceding 
case already developed in (33). Thus, for any triple mirror, the 
order of reflections being I, 2, 3, the resultant operator ft *=ft 123 is 



. n 2 + 31 n 3 . nj - 8a ia a 23 n 3 . n^ (37)* 

Here again, the second term being symmetrical, but not the third, 
nor the last, I2 132 , etc., will differ from 12 12D . In other words, an 
incident beam of parallel rays r will give rise to six reflected beams 
represented by 

etc., 



which will, generally speaking, deviate from one another.- These 
six reflected beams become parallel to one another, i.e. 12 becomes 
independent of the order of the three successive reflections when 
and only when 

U !SS *!8 rB *81 S!BO 

abolishing the third and the fourth terms in (37), i.e. when the 
three component mirrors are mutually perpendicular. Moreover, 
n i> n 2> n s being in that case a triad of normal unit vectors, the 
dyadic n a . nj+n 2 . n 2 +n 3 . n 3 becomes an idemfactor or I, and 
the operator is reduced to 

12= -i ; /. r'=-r. (37*) 

That is to say, every incident ray, independently of the order of its 
three reflections, is sent back parallel to itself. The orthogonal 
triple mirror is a central mirror. Such triple mirrors, in fact, are 
used in practice, their orthogonality being now attainable to 
a high degree of approximation. 

Returning to the general triple mirror, write r,, r 2t r z for the 
direction cosines of the incident ray or its projections upon the 
three normals n K . Then the direction cosines of the reflected 
ray r' will be linear homogeneous functions of r,, r 2 , r 3 , viz., by (37), 

*i - [ - 1 + 4V + V - S^^aiK ~ Ni2 + 4^a 13 ]r 2 - 2a 13 r 3 , (38) 
and so on. 

* This may, at first sight, seem very complicate. Try, however, to treat the* 
general triple mirror by the ordinary methods. Then (37) will be found very 
simple. 
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Exercises. Write down, by the general formula (37), the remaining 
two direction cosines r,', rj, thus obtaining three scalar formulae 
with n lf n a , n a as reference system. 

Assuming that n lf n s , n 3 are not coplanar, i.e. that the three 
mirrors constitute a proper pyramid (not a prism), any incident ray r 
can be put into the form 



where p lt etc., are some scalars. Determine them in terms of the 
orthogonal projections ^=1^, etc. (Notice that pi=r v etc., when 
and only when the mirrors are mutually perpendicular.) 
Similarly determine the coefficients of 



and replace the formulae (38) by three equations giving the p K ' in 
terms of the p K , and in terms of the properties of the mirror, i.e. a M , 
etc. Discuss the general formulae thus obtained, and apply them to 
a number of particular cases. 

Treat the prismatic arrangement of the three component mirrors, 
when HI, n a , n 3 are coplanar ; i.e. such that an x +6n a +cn 8 =o, a, 6, c 
being three scalars, and so on. Consider first an incident ray JL to 
the edges (all represen table by one e), and then any oblique or skew 
ray. In the former case imagine the mirror pierced at a place, just 
to let through the incident ray ; or else let the mirrors consist of 
non-silvered (plane parallel) glass plates. 

The angle ^ between the incident and the finally reflected ray 
is given by cos^=rJ2r, i.e., according to (37), by 

J(i - cos $) = (rf + r + r 3 2 ) - 2[a 12 r 1 r 2 + etc.] +.4iaVi- (39) 
Each of the two bracketed expressions is symmetrical with respect 
to the suffixes I, 2, 3, but the last term is not. Thus, if we replace 
the order of reflections 123, by 132, say, and call ^' the angle 
between the corresponding finally reflected ray and r, we have the 
elegant formula 

cos ^' - cos ^ = 8r 1 a 23 (a 12 r 3 - a 13 r 2 ) . (40) 

Similarly for any other permutation. Notice that 321 gives the 
same reflected ray as 1^3. 

If the three mirrors constitute a regular pyramid^ i,e. if 

M = *-- cos , say, 
then formula (39) becomes 

}( I - cos $ - (rf + r 2 2 + r 3 2 ) - 2 cos eofor, + *y 3 + r^) 

. (390) 
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More especially, if the incident ray r is equally inclined to the three 
reflecting planes of the regular pyramid, we have 



and therefore, remembering that I - cos ^ = 2 sin 2 -, 

sin-=cot-/v/i -2 cosw+- cos 2 w. (396) 

In this case the angle \f/ is, of course, independent of the order of 
reflections. The reflected beams (corresponding to a beam of 
parallel incident rays) are not parallel to one another, but are 
equally inclined to and symmetrically disposed around the direction 
of the incident beam. These reflected beams will coincide in 
direction when and only when <o = 90, i.e. when the mirror becomes 
an orthogonal and, therefore, a central mirror. In fact, for 

ci>=-, formula (39^) gives sin- = i, i.e. ^ = 180, which means a 

reversal of the incident ray. This, as we already know, is the 
property of an orthogonal mirror with respect to any incident ray. 
Quadruple and more Complicated Mirrors. Multiplying (37) by 
I - 2n 4 . n 4 , where n 4 is the normal of the fourth mirror, the reader 
will obtain the operator ft=12 W34 , and similarly for quintuple 
reflections, and so on. As was already remarked, some of the 
normals n K may reappear, that is, we may have, for instance, a 
quintuple reflection from but two mirrors ; if such be the case, 
it is enough to put n 3 =n 6 =n 1 , and n 4 =n 2 . 

As an example of this kind the reader may treat the rotation and 
progress of a skew ray incident upon an interior face of an equilateral 
prism built up of three plane mirrors ; n lf n a , n$ coplanar, n 4 =ni, 
n a =n 2 , and so on. 

Reversal of the Order of Reflections. For any multiple mirror, let 
r be the incident ray, r' the finally reflected ray when the order of 
reflections is 123 ... N - 1, N, and T" the finally reflected ray when 
the order of reflections is reversed, i.e. N, N- i, ... 321. Then, 
ft being ft*fty-i ... ^y^^it as on P a S e 2 5> ^ e h a ve 

r'=ftr and r"=rfl, 
whence 

rr ' = rfir, r "r = rr " = rftr, 

and therefore, for any multiple mirror, 

r'rr"r, (41) 
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while rV^rft^r. That is to say, the reflected rays r', r", although 
not parallel to one another, are always equally inclined to the incident 
beam r. The equality ^123 = ^321 mentioned on page 30 is but 
a special instance of this general property. 



8. Dyadic representing the most General Reflector. 

Consider the successive reflections of a ray r from any number 
of mirrors, whose normals are n x , n 2 , etc., and let ft, with the 
order of reflections subintended, be the corresponding resultant 
operator or the " reflector " belonging to the whole, ordered, 
system. (This comprises also the case of any number of curved 
mirrors ; the only difference being that when they are plane, all 
the normals are immediately given, and when curved, the successive 
normals are to be found by the application of the " transfer for- 
mula," developed in a previous Section.) Now ft being the product 
of dyadics of the form 

i - 2n . n, 

is itself a dyadic. And since our space has but three dimensions, 
all the consequents * of any reflector, no matter how complicated 
the system to which it belongs, can be expressed as linear homo- 
geneous functions of three non-coplanar vectors, which can always 
be made unit vectors, let us say, f, g, h. The antecedents belonging 
to these consequents will then, for any given multiple mirror, be 
certain three vectors, let us say, a, b, c. Thus, the most general 
reflector will have the form of a trinomial dyadic, 

ft =a . f +b . g +c . h. (42) 

The plain meaning of this is that any incident ray r will be reflected, 
in the prescribed order, as 

r' = a(fr) +b(gr) +c(hr). (420) 

The consequents f, g, h need not necessarily constitute a normal 
system ; it may sometimes be convenient to choose an oblique 
system of vectors. The choice will depend upon the properties of 
the system of mirrors. 

The above result may also be stated by saying that the most 
general ft is a linear vector operator, or r' a linear homogeneous 
vector function of r, as was already mentioned. 

* According to Gibbs' terminology, the first vector, a, in a dyad a . b is called the 
antecedent, and the second vector, b, the consequent. 
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Splitting into Cartesians, for instance, and counting up the 
coefficients expressing the components of r' in terms of those of r, 
the reader will see at once that the most general linear vector 
operator implies nine independent scalar or numerical data. Let 
us count the number of such data, implied in the most general 
reflector 12, upon the formula (42) itself. 

The consequents, f, g, h, which can always be taken as unit 
vectors, and, moreover, always as a normal system of unit vectors,* 
imply 2 H- 1+0= 3 scalar data (for the direction f being fixed, g 
requires but one datum, and h is codetermined, apart from the 
irrelevant sign or sense). Now it is always (i.e. for any trinomial 
dyadic) possible to give to the normal system f, g, h such an orienta- 
tion as to convert the antecedents also into some normal system of 
vectors, t but, of course, not unit vectors, after f, etc., have been 
made so. Thus the orientation of the system of antecedents 
implies three more data, and the size or absolute value of each 
antecedent one datum. Thus the total number of scalar data 
implied in the most general trinomial dyadic is nine. 

And since our reflector 12 is such that it leaves the size of any 

operand unaltered, i.e. 

r' 2 =r 

for any r, the number of independent scalar coefficients contained 
in the most general reflector 12 is reduced to eight. These coefficients 
depend upon the properties of the system of mirrors and, of course, 
of the prescribed order of successive reflections. No matter, 
however, what their arrangement and how many of simple mirrors 
are built up into the system, all of its properties are crowded into 
eight coefficients, corresponding to each prescribed order. On 
permutating the order, the operator, or its coefficients, will, in 
general, be changed ; in some special cases they may remain 
intact. The reversal of order has, however, the perfectly general 
property expressed by (41). 

A dyadic, or linear vector operator, such that J2r=ri2, for any r, 
is called a self-conjugate or symmetrical one. ^ The simple reflector 

* Although this may sometimes be not the most convenient choice. This, how- 
ever, does not influence the result of our counting. 

t A simple proof is given in J. W. Gibbs* lectures on Vector Analysis, edited by 
E. B. Wilson, New York, 1907, p. 304. Our dot, a separator, is not to be con- 
founded with Gibbs' dot, which is the symbol of scalar multiplication. The 
notation adopted in the present volume is essentially that of Oliver Heaviside. 
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i - 2n . n is, obviously, such an operator. But that corresponding 
to the general double mirror is no longer self-conjugate. Other 
details of the discussion of the multiple reflection formulae may 
be left to the reader. As to further examples and exercises, 
besides those proposed, the reader will doubtless invent them himself 
without difficulty. 

9. Differential Properties of Refracted or Reflected Pencils. 

Although the chief object of this book, viz. the tracing of a single 
ray, is in its essence exhausted by what has been said in the pre- 
ceding sections, it may be well to add one more section in order 
to give the reader at least some cursory hints about the treatment 
of differential properties of beams or pencils of rays sent through an 
optical system : properties, that is to say, which concern the 
passage from one ray of a pencil to its neighbours, or the passage 
from a given wave-length, and the corresponding index /*, to an 
infmitesimally different one (astigmatism, chromatic aberration). 

Problems of this kind do not, as a matter of fact, require anything 
essentially different from what has already been given in the 
earlier part of the book. They are solved, in each case, by simply 
differentiating both sides of the vector formulae concerning a 
single ray, obtainable by the above method. Now, the differentia- 
tion of vector expressions does not differ from that of scalar ones, 
provided of course that we differentiate with respect to a 
scalar magnitude as the independent variable or parameter. As 
a precaution hardly anything more is needed than the preservation 
of order of the (vector) factors of a vector product and, in case 
of dyadics, the discrimination between antecedent and consequent. 

Thus, if f means the " differential coefficient " or the flux of r 

9r 
with respect to any scalar variable x, i.e. if r stands for ^-, the 

flux of a sum a+b of two vectors is a+b, and that of the scalar 
product ab is ab +ab. Similarly, for the vector product, 



and so on. Or, using " differentials, 1 * 

d(a + b) = rfa + db, d(ab) 
where rfa, db are infinitesimal vectors. Similarly 
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> 

the second term being negatived to compensate for the reversal of 
the order of its two factors. Further explanations, if at all needed, 
will be found in any treatise on Vector Analysis.* 

Such being the case, it will be enough briefly to consider the 
passage of a pencil of rays r through a single refracting surface 
(reflection being, formally, but a sub-case). Thus, returning to 
the fundamental refraction formula (II), Section I, we have for 
any ray of an incident pencil, let us say for the " central " ray of 

the pencil, 

s' = s + n, (II) 

where n is the normal of the refracting surface, of any shape (free 
of singularities), at the point struck by the central incident ray s. 
Let the refractive indices /*, \L be fixed, and let us pass from the 
central ray s (to which corresponds the central ray s' of the refracted 
pencil) to an infinitesimally neighbouring ray s-f ds. If all the 
rays of the pencil are parallel, then, of course, ds=o, s representing 
the whole pencil. But in general we shall have d&o. At any 
rate, since s 2 =/x 2 , the relation srfs=0, or rfsXs, will hold. The 
corresponding refracted ray will be s'-fds', and da' Is' as before. 
Now, passing from one ray to another means, in general, varying the 
normal n, and therefore also g. Thus, subtracting the original or 
central s' from the new vector s'+rfs', we have for the required 
change of the refracted ray, 

da' =? dz + g dn + n dg, (43) 

where, by the definition of g, p. 2, 

dg=p sin i . di - // sin i'. di' t 
and by the sine-law, /*' cos i'. di' =/* cos i . di, so that ultimately 



f being a known function of i alone. Notice that, n being a unit 
vector, we have always ndn=o. In the case of reflection it is 
enough to replace // by -/*. 

In essence, formula (43) already contains the -answer to all 
thinkable questions concerning the properties of the infinitesimal 
refracted pencil in terms of those of the incident one. It is enough 
to read it intelligently and to develop it appropriately in each 
particular case. 

* Cf., for instance, pp. 22-48 of my Vectorial Mechanics. Only a small part of 
what is there given will be required for the purposes of geometrical optics. 
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Let , v be any curvilinear coordinates defining the position of 
a point at the refracting surface o-. Then " incident pencil given " 
means da given as a function of w, v throughout an infinitesimal 
region round U Q , v , the point struck by the central ray. Now, the 
surface <r being given (round that place at least), say, through the 
expression of its squared line-element, 

dP = E du 2 + 2Fdudv + G dv\ (44) 

dn is also known as a function of the two coordinates. The same 
is true of dg, which differs from di by a known factor f(i). 

Thus (43) finds ds' as function of u, v round w , V Q ; that is, it 
gives all the properties of the refracted pencil. For s', n or s', 
w o> V Q gi ye t ne direction and the position of the central pencil, and 

s'+rfs' with Uo+du, V + dv 

localize and determine the direction of any ray of the refracted 
pencil. Thus every question concerning the latter pencil is reduced 
to purely geometrical constructions. 

Employing u, v (for instance, <f>, for a sphere) as our independent 
variables, we have, of course, 



and in quite the same way 



If the incident pencil consists, for instance, of parallel rays, so 
that ds=o, then 

(45) 



whence nda'=*dg, and Vn^s'=^Vnrfn, two convenient formulae 
replacing (45). The chief question in any such investigation is 
whether the neighbour ray passing through P(u Q + du, v Q + dv) 
meets at all, at finite or infinite distance, the central refracted ray 
passing through P (u , V Q ). In order to put this question in 
vectorial language, ifttroduce the infinitesimal vector dl, drawn on 
the surface from P to P, its size dl being determined by the quad- 
ratic form (44) with E, F, G as given functions of u, v. Then the 
said rays meet or not (are coplanar or not) according as the scalar 
magnitude 
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vanishes or not. Remembering that Vfl's'=0, identically, this 
criterion-expression is reduced to 



which, by (II) and (45), is 

Q [dg Van +gVs dn +g 2 Vn dn]dl. (46) 

Given the refracting surface and the incident pencil s, together 
with the central normal n, this expression will vanish or not accord- 
ing to the choice of the point P, i.e. of the vectorial line element ifl. 
In general we shall have j2=o only for certain particular directions 
of dl upon the interface. Having established the conditions 
for (2=0, it will be a comparatively easy task to localize the point 
of intersection of the P-ray with the central ray of the refracted 
pencil, and thus to build up the caustics. If the incident beam 
is normal, then, whatever the shape of the refracting surface, 
the most convenient curvilinear coordinates u t v will be those 
measured along the lines of curvature of that surface drawn through 
P as origin, in which case (44) is reduced to dP=*Edu* + Gdv 2 . 
For oblique incidence, however, another network may in certain 
cases be found more convenient. 

As an exercise the reader may treat the case of normal incidence 
upon any refracting surface, such as an ellipsoid, a cylinder or 
a cone ; and the case of oblique incidence for a spherical surface, 
thus deriving vectorially the well-known formulae for the sagittal 
and the meridional section of a pencil. 

Since these and similar problems have nothing to do with " ray- 
tracing, M the proper argument of this little volume, we cannot 
enter upon them any further. The last mentioned, and similar, 
exercises may, however, be warmly recommended to the reader 
as good and interesting opportunities of training in the use of 
vector language. 
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